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Abstract. An intrinsic curvature model is investigated using the canonical Monte 
Carlo simulations on dynamically triangulated spherical surfaces of size upto N — 4842 
with two fixed-vertices separated by the distance 2L. We found a first-order transition 
at finite curvature coefficient a, and moreover that the order of the transition remains 
unchanged even when L is enlarged such that the surfaces become sufficiently oblong. 
This is in sharp contrast to the known results of the same model on tethered surfaces, 
where the transition weakens to a second-order one as L is increased. The phase 
transition of the model in this paper separates the smooth phase from the crumpled 
phase. The surfaces become string-like between two point-boundaries in the crumpled 
phase. On the contrary, we can see a spherical lump on the oblong surfaces in the 
smooth phase. The string tension was calculated and was found to have a jump at 
the transition point. The value of a is independent of L in the smooth phase, while 
it increases with increasing L in the crumpled phase. This behavior of a is consistent 
with the observed scaling relation a ~ (2L/N) U , where v ~ in the smooth phase, 
and v = 0.93 ± 0.14 in the crumpled phase. We should note that a possibility of a 
continuous transition is not completely eliminated. 
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1. Introduction 

Phase transition of surfaces has been investigated to understand string models and 
biological membranes [H [21 El HI El [6], [7] . Efforts has been devoted to clarify the phase 
structure [El EJ [101 [111 [12] m the model of Helfrich, Polyakov, and Kleinert (HPK), which 
is defined by a two-dimensional curvature Hamiltonian [131 El 03] • The self-avoiding 
property of surfaces are crucial for biological membranes [161 El EEE], however, it is 
less important for string models pQ. A surface model that is allowed to self- intersect is 
conventionally called a phantom surface model, which has long been studied numerically 
on tethered surfaces p21 EH El] and on fluid surfaces (^ESIESESEEIETIEEIEHIEO]. 
Recent numerical simulations on phantom HPK models revealed that there is a first- 
order transition between the smooth phase and the crumpled phase [3H [321 E31 El] . The 
results in [331 [31] imply a possible phase transition in biological or artificial membranes, 
because the Hausdorff dimension H is H < 3 in the crumpled phase even though the 
surface is allowed to self- intersect. 

The string tension a is a key observable to understand the phase structure of 
surfaces [24]. If a surface becomes oblong, string-like, and of length 2L, the free energy 
of the surface can be expressed by a2L. Therefore, a can be extracted from the grand 
canonical ensembles of surfaces with fixed boundaries of sufficiently large distance 2L. 
It was reported that an extrinsic curvature model undergoes a first-order transition, 
where the string tension becomes non-zero in the smooth phase close to the transition 
point in the limit of L —>■ oo [28, 29J. The results of the extrinsic curvature model in 
[281 129] show that the order of the transition changes depending on L. 

However, it is still unknown how the distance 2L between the boundaries influences 
the phase transition of an intrinsic curvature model on dynamically triangulated 
surfaces. The intrinsic curvature model was first studied by Baillie et.al. [35l [36l [371 [38] . 
A surface model with deficit angle term, which is an intrinsic curvature model, has a first- 
order transition on tethered surfaces [39j . Another intrinsic curvature model undergoes 
a first-order transition, which is independent of the surface topology [40j HH [42] . It was 
also reported that the order of transition changes from first to second as L increases from 
L = L to L = 2L , where L is the radius of the original sphere constructed to satisfy the 
relation S\/N = 1.5 [13]. However, the string tension was unable to be extracted from 
the simulations on the tethered surfaces, because the surfaces do not always become 
string-like even when L is relatively large compared to L$. On the contrary, the fluid 
surface can be oblong in the limit of L — > oo. Therefore, we expect that not only the 
dependence of the transition on the distance 2L but also the string tension a is evaluated 
from the intrinsic curvature models on fluid surfaces with the fixed boundary-vertices. 

In this paper, we will study the dependence of the transition on the distance 
2L of the boundary of spherical fluid surfaces by MC simulations with dynamical 
triangulations. The string tension can be obtained by the canonical MC technique, 
although it is originally defined by the grand canonical ensemble. We confine ourselves 
into a phantom surface model in this paper. However, it is natural to expect the surfaces 
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are not completely self-intersected even in the crumpled phase, because they span a long 
distance between two point-boundaries and consequently become stretched out. 

2. Model 

The Gaussian bond potential Si and the intrinsic curvature energy S3 are defined by 
S 1 = ^(X J -X J ) 2 , S 3 = -£log(*/27r), (1) 

(ij) i 

where i n Si is the sum over bonds (ij) connecting the vertices i and j. The bonds 
(ij) are the edges of the triangles. Si in Eq. (CD) is the vertex angle, which is the sum 
of the angles meeting at the vertex i. We call S3 the deficit angle term, because 5— 2tt 
is just the deficit angle. We note that J2i(o~i~ 2vr) is constant on surfaces of fixed genus 
because of the Gauss-Bonnet theorem. 

The partition function of the fluid surface model is defined by 

, N 

Z = H \{dx l eM-s(x 1 r)}, (2) 

r J i=i 
S(X, T) = Si + «S 3 , 

where a is the curvature coefficient, and YjT denotes the summation over all possible 
triangulation T. /' in Eq.(j2J) denotes the boundary condition in which two vertices 
are fixed and separated by the distance 2L(N), which will be discussed in the following 
section in more detail. S(X, T) denotes that the Hamiltonian S depends on the position 
variables X of vertices and the triangulation T. 

We should comment on the unit of physical quantities in the model. We can express 
all physical quantities by unit of length in terms of a, which is a length unit in the model. 
Hence, the unit of Si and Lq(N) are a 2 and a, respectively. We fix the value of a to 
a = 1 in this paper, because the length unit can be arbitrarily chosen on the basis of the 
scale invariant property in the model. The unit of a is expressed by kT, where k is the 
Boltzmann constant, T is the temperature. Note also that varying the temperature T 
is effectively identical with varying the curvature coefficient a in the model. 

We note also a relation between the deficit angle term S3 and the integration 
measure HidXiq^ 2 in the partition function of a fluid surface model [H], where % 

3 /2 

is the coordination number of the vertex i. The weight factor q i ' can also be written 
as exp[(3/2) J2i l°g9t]- Then, replacing by 5i and 3/2 by a the continuous number, 
we have exp[a Y^i log5j]. Introducing the constant term into the Hamiltonian, we have 
S 3 in Eq.©. 

3. Triangulated spheres and Monte Carlo technique 

Triangulated surfaces are obtained by dividing the icosahedron. By splitting the edges 
of the icosahedron into ^-pieces, we have a surface of size N = 10£ 2 + 2. These surfaces 
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are characterized by N$ = 12 and N 6 = N — 12, where N q is the total number of vertices 
with a co-ordination number q. 

The spherical surfaces are constructed by fixing the radius Lq(N) so that the 
Gaussian potential S\/N is approximately equal to Si/N~3/2. We assume the following 
values of L (N): L (N) = 7 for N= 1442, L (N) = 9.3 for iV = 2562, L (JV) = 12.8 for 
iV = 4842. Under those values of L (N), the relation Si/N— 1.5 is almost satisfied on 
those spheres. Note that the value of the radius L (N), which satisfies S±/N~3/2, may 
slightly change depending on how to construct the sphere. 

The distance 2L(N) between the boundary vertices is fixed to three different values 
such that 

L(N) = 2L (N), 3L (N), 5L (N) (3) 

where L (N) is the radius of sphere constructed from the icosahedron as described 
above. 



(a) L(N)=2L (N) (b) L(N) = 3L (N) (c) L(N) =5L (N) 

iV = 4842 iV = 4842 iV = 4842 




Figure 1. Starting configurations under the conditions (a) L(N) = 2Lq(N), (b) 
L(N) = 3L Q (N), and (c) L(N) = 5L (N). The symbols + at the terminal points 
represent the position of the boundary vertices. The distance between the boundaries 
is 2L(N) in each surface. These figures are drawn in scales which are different from 
each other. 



Figures [T^a), [T^b), and[]Tc) show the starting configurations of the MC simulations 
for L(N) = 2L (N), L(N) = 3L (N), and L(N) = 5L (N). In order to show the position 
of boundary points, the symbols + are plotted at the points in each figure. These figures 
are drawn in scales which are different from each other. The size of spherical lump at 
the center of each snapshot is a little bit smaller than that of the original sphere of 
radius L (N). 

The surfaces in FigsJT](a)-[T](c) are obtained by pulling the boundary vertices out of 
the original spheres. We choose the boundary points and fix the positions as follows: 
Firstly, we choose a pair of vertices, which are on a straight line passing through the 
center of the sphere. The canonical x-coordinate axis in R 3 = {(x, y, z)\x, y, z G R} is 
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chosen along the straight line. The distance between these two vertices on the axis is 
given by 2L (N) at the beginning. Secondly, the distance is enlarged from 2L (N) to 
2L(N) along x-axis during the thermalization MCS (Monte Carlo Sweep), where L(N) is 
given in Eq. ([3]). The vertices are pulled away from the sphere during the first 0.25xl0 7 
thermalization MCS, and then 0.75xl0 7 or more thermalization MCS is performed after 
the expansion was finished. Equilibrium statistical mechanical conditions seem to be 
violated slightly in the first 0.25 xlO 7 MCS because of this forced expansion of vertices. 
However, it should be noted that these MCS are iterated only for the constructions 
of starting configurations. These MCS are neither the production MCS for sampling 
nor the thermalization MCS. Therefore, the final results are independent of the method 
to get the starting configurations because the number of thermalization MCS is large 
enough. 

The distance 2L(N) is obtained by enlarging the distance between the boundary 
points from 2L(N) = 14 to 2L(N) = 28 on the N = 1442 surface for example. Two 
fixed vertices are shifted by distance 7 x 10~ 5 at every 25 MCS to the opposite direction 
to each other along x-axis during the first 0.25 x 10 7 MCS. We must emphasize that 
the results of the simulation is completely independent of how the surface with the 
boundaries of distance 2L(N) is constructed from the sphere of diameter 2L Q (N). 

The distance between boundary vertices on the surface of L(N) = 5L (N) is just 
10Lo(iV) for example. The boundary conditions are thus imposed on the model by these 
two fixed vertices of distance 2L(N). If it were not for the boundary condition, we have 
S\/N = 3/2, which comes from the scale invariant property of Z. However, as we will 
see later, this relation is slightly broken due to the boundary condition. 

The vertices X are shifted so that X' = X + 5X in MC, where 5X is chosen 
randomly in a small sphere. The new position X' is accepted with the probability 
Min[l, exp(— AS)], where AS = S(new) — 5 (old). We use a sequence of random numbers 
called Mersenne Twister [45J for three dimensional random shift of X and for the 
Metropolis accept /reject. The radius of the small sphere for the shift 6X is chosen 
so that the rate of acceptance for X is about 50%. We introduce the lower bound 
1 x 10 s for the area of triangles. No lower bound is imposed on the bond length. 

The dynamical triangulation is performed by the standard bond flip technique 
for integrating the degrees of freedom T in the partition function Z of Eq.(j2J). The 
bonds are labeled with sequential numbers, where the total number of bonds Nb is 
given by Nb = 3N — 6. The odd-numbered bonds are firstly chosen to be flipped, 
and secondly the remaining even-numbered bonds are chosen. The flip is accepted 
with the probability Min[l, exp(— AS)]. iV-updates for X and Nb/2- updates for T are 
consecutively performed and make one MCS. 

In order to improve computational efficiency, we renumber the vertices in order close 
to one boundary vertex at every 0.5 x 10 5 ~ 5 x 10 5 MCS. Vertices freely diffuse over the 
surface because of the fluidity introduced by the bond flip in dynamically triangulated 
MC simulations. As a consequence, the vertex number is eventually distributed at 
random over the surface. The computations are very time consuming on such random 
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surfaces. This is typical of the computations on dynamically triangulated fluid surfaces. 
However, the renumbering the vertices can save the computer time. The speed of 
computation is increased about 100% or more by the renumbering on the N = 4842 
surface. 

4. Results 

Monte Carlo simulations were performed on the surfaces of size iV = 1442, iV = 2562, and 
iV = 4842. The distance between the boundaries is 2L(N)=4L (N), 2L(N) = QL (N) } 
and 2L(N) = 10L (N) on each surface. The total number of MCS is about 0.8 x 10 8 ~ 
1.5 x 10 8 at transition point a c after the thermalization MCS. Relatively small numbers 
of sweeps were performed at non-transition points a=/a c . 



(a) a = 1850, L(N)=2L (N) (b) a = 1900, L(N) = 2L (N) 




Figure 2. Snapshots of the L(N) — 2Lq(N) surface obtained at (a) a — 1850 (crumpled 
phase) and (b) a = 1900 (smooth phase). The symbols + denote the position of the 
boundary points. The snapshots are drawn in the same scale. The surface size is 
iV = 4842. 

Firstly, we show snapshots of the L(N) = 2L (N) surfaces obtained at a = 1850 
and a = 1900 in FigsJ2(a) and 12(b), respectively. The distance between the boundaries 
is 2L(N) = 4Lq(N). The surface size is N = 4842 in both snapshots. The surface in 
FigJ2(a) is crumpled and is typical of the crumpled phase, while the surface in FigJ21(b) 
is partly swollen and is typical of the smooth phase. The crumpled surface seems linear 
rather than branched-polymer, while the smooth surface has a lump of smooth sphere, 
which stays outside the boundary points and touches one boundary point. Almost all 
vertices are included in the lump. We confirmed that the lump is a smooth surface by 
slicing it into small sections. 

Figures [3(a) and [3(b) are the snapshots obtained at a = 1900 and a = 1950 with the 
distance 2L(N) = 6L (N), where iV = 4842. Both figures are drawn in the same scale. 
The snapshot in Figj3(a) is the one in the crumpled phase, while that of Figj3(b) is the 
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(a) a = 1900, L(N)=3L (N) 0) a = 1950, L(N) = 3L Q (N) 




Figure 3. Snapshots of the L(N) — 3Lo(N) surface obtained at (a) a — 1900 (crumpled 
phase) and (b) a= 1950 (smooth phase). The symbols + denote the position of the 
boundary points. The snapshots are drawn in the same scale. The surface size is 
iV = 4842. 

one in the smooth phase. We can also see that the crumpled surface is almost linear 
and the smooth surface includes a smooth lump as in Fig. [2(b) for L(N) = 2L Q (N). The 
lump in the smooth surface is located outside the boundary points just as in Fig. [2](b) 
for L(N)=2L (N). 



(a) a = 2000, L(N)=5L (N) (b) a = 2050, L(N) = 5L Q (N) 




Figure 4. Snapshots of the L(N) = 5Lq(N) surface obtained at (a) a — 2000 (crumpled 
phase) and (b) a = 2050 (smooth phase). The symbols + denote the position of the 
boundary points. The snapshots are drawn in the same scale. The surface size is 
iV = 4842. 

Snapshots obtained at a = 2000 and a = 2050 are shown in Figsfjja) and BJb), 
where the distance between boundaries is 2L(N) = 10L (N), and the surface size is 
iV = 4842. The crumpled surface in FigJUa) is apparently linear, and the smooth 
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surface in FigJH(b) includes a spherical lump as those in the cases L(N) = 2L (N) and 
L(N) = 3L (N) in Figj2](b) and Figj3^b). Contrary to the previous cases in FigsJS^b) and 
&h), the spherical lump stays between the fixed boundaries in FigQJb). However, the 
lump is not yet trapped at either boundary. Surfaces span inside the fixed boundaries 
in both crumpled and smooth phases in the case L(N) =5Lq(N). 
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Figure 5. The bending energy S^/Wb vs. a obtained under the conditions (a) 
L(N) = 2L (N), (b) L(N) = 3L (N), and (c) L(N)=5L (N). N B is the total number 
of bonds. Solid lines were drawn by the multihistogram reweighting technique. 



The bending energy S2 is defined by using the unit normal vector rij of the triangle 
i such that 

5 2 = ^(l-ii,.^), (4) 

where YUj denotes the sum over all nearest neighbor triangles i and j that share the 
common bond. S2 is not included in the Hamiltonian, however, it reflects how smooth 
the surface is. 

Figures Eta)-[5tc) show S2/NB against a obtained on the surfaces of L(N) = 2L (N), 
L(N) = 3L (N), and L(N) = 5L (N). Nb is the total number of bonds. Solid lines were 
drawn by the multihistogram reweighting technique [415] . The reweighting is done by 
using 4 ~ 5 data close to the transition point. The dashed lines, which simply connect 
the data, are drawn to guide the eyes. We can obviously see that S 2 /N B disco ntinuously 
changes against a in all three conditions for L(N). A phase transition can be called 
a first-order one if some physical quantity discontinuously changes. Therefore, the 
discontinuous change of S2/NB in Fig[5|a)-[5](c) indicates that the model undergoes 
a first-order transition under the conditions L(N) = 2Lq(N), L(N) = 3Lq(N), and 
L(N) = 5L (N). 

The standard deviations in S^/iV^ are calculated and plotted in the figures as error 
bars. However, the errors are small and hardly to be seen. In the following, the standard 
deviations are also calculated and plotted as error bars in all the figures just as those in 
FigO 
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Figure 6. The mean square size X 2 vs. a obtained under the conditions (a) 
L(N) = 2L (N) , (b) L(N) = 3L (N) , and (c) L(N) = 5L (N). Dashed lines drawn 
vertically denote the phase boundary between the smooth phase and the crumpled 
phase. 



The mean square size X 2 is denned by 

x 2 = ^£(^-*) 2 > x=4£*<, (5) 

i i 

where X is the center of the surface. X 2 reflects the size of surfaces even on surfaces 
with fixed boundaries of distance L(N) given by Eq. (j3J). In the crumpled phase, 
X 2 is expected to be almost constant against a, because the surface spans from one 
boundary point to the other and is almost linear, and therefore there is no reason for 
large fluctuations of X 2 . On the contrary, X 2 largely fluctuates in the smooth phase 
because it depends on the position of the spherical lump, which can be seen in the 
snapshots shown in Figs. EJ^b), E^b), and|H(b). X 2 is expected to remain small when 
the lump stays in the middle of the two boundary points, while it becomes relatively 
large if the lump moves apart from the center of the boundary points. 

Figures ma)-E](c) show X 2 against a for L(N) = 2L (N), L(N) = 3L {N), and 
L(N) = 5L (N). The dashed lines drawn vertically in the figures denote the phase 
boundary between the crumpled phase and the smooth phase. We find from the figures 
that an expected large fluctuation can be seen in X 2 of the N = 4842 surfaces at the 
smooth phase. Whereas X 2 of the N= 1442 surfaces remains almost constant against a. 
The reason of this is the size effect. The total number of vertices in the spherical lump, 
which appear on the surface in the smooth phase, is relatively small, and therefore X 2 
is not influenced by where the lump is in such small sized surfaces. 

We plot the Gaussian bond potential Si/N against a in FigsJ7ta)-[7tc). Si/N 
should be equal to 3/2 whenever no specific boundary condition is imposed on the 
model as mentioned previously. However, Si/N in Figsd(a)-0(c) obviously deviates 
from 3/2. The reason of this is because the surface is fixed with two vertices separated 
by the distance 2L(N). The figures were plotted in the same scale along the vertical 
axis. We find from the figures that the deviation of S\/N from 3/2 increases with 
increasing L(N). It is also clear that S±/N discontinuously changes in FigsJ7](a)-[7](c). 
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Figure 7. The Gaussian bond potential Si vs. a obtained under the conditions (a) 
L(N) = 2L {N) , (b) L(N)=3L (N) , and (c) L(N) = 5L (N). Solid lines were drawn 
by the multihistogram reweighting technique. 



The discontinuity of S\/N indicates that the model undergoes a first-order transition 
under L(N) = 2L Q (N), L(N) = 3L (N), and L(N) = 5L Q (N), in accordance with the 
prediction from the behavior of S^/iVe in Figsj5](a)-[5](c). 

The shape of surfaces is linear and seems string-like when the surfaces are in 
the crumpled phase as we have seen in the snapshots. Moreover, the surfaces are 
considered to be sufficiently oblong even when they are in the smooth phase except 
under L(N) = 2L Q (N). In fact, the size of the spherical lump in the smooth phase for 
the L(N) = 2Lq(N) surface in FigJ2](b) is comparable to the half distance between the 
fixed boundaries, and therefore we are unable to see the surface as a string-like one. On 
the contrary, the lump of sphere appears sufficiently small compared with the distance 
2L(N) on the surfaces of L(N) = 3L (N) and L(N) = 5L (N) in FigsJ^b) andHfb), 
and therefore it is possible to see the surfaces as string-like when L(N) = 3L (N) and 

L(iV) = 5Lo(iV). 

Therefore, we expect [21] 

Z(a;L) ~exp(-(j2L) (6) 

when L(N) = 3Lq(N) and L(N) = 5Lo(N). Then, by using the scale invariance of the 
partition function, we have [H [21] 
2St - 3N 



a 



2L 



(7) 



where N and L remain constant throughout the simulations. From this expression, we 
simply expect that a discontinuously changes at the transition point, because S\ has a 
jump as we have seen in FigsJ7](a)-[7[c). 

Figures E](a)-[H](c) are plots of a against a under the conditions L(N) = 2L Q (N), 
L(N) = 3L (N), and L(N) = 5L (N). Although a is not well-defined under L(N) = 
2Lq(N) at least in the smooth phase because the surface is not always string-like as 
can be seen in the snapshot in Fig 121(a), it is plotted in FigJEIa). On the contrary, a is 
considered to be well-defined in FigsJH](b) and[8](c), and we see that o has an expected 
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Figure 8. The string tension a vs. a obtained under the conditions (a) L(N) — 
2L (N), (b) L(N) = 3L (N), and (c) L{N) =5L {N). Solid lines were drawn by the 
multihistogram reweighting technique. 



behavior against a under L(N) = 3Lo(N) and L(N) = 5Lq(N); a discontinuously changes 
at the first-order transition point. 

We also find from the results of the N = 4842 surfaces shown in FigsJH](b) and W^c) 
that: 

(i) a remains almost constant against a in the smooth phase 

(ii) a is independent of L(N) in the smooth phase. 

(iii) a increases with increasing L(N) in the crumpled phase 

On the N = 1442 surfaces in FigsJHl(b) and[8](c), we find that a increases with increasing 
iV in the smooth phase. However, we consider this is only due to the size effect. 
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Figure 9. Log-log plot of a against 2L/N, obtained at the smooth phase close to the 
transition point and at the crumpled phase close to the transition point. The straight 
lines are drawn by fitting the data to the relation of Eq.©. 



In order to see the scaling behavior of a [21] 
docf, t = 2L/N 



(8) 
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we plot cr versus t = 2L/N in FigJH] in a log- log scale, where v in Eq.flS]) is some 
scaling exponent. The empty (solid) symbols denote that a is obtained in the smooth 
(crumpled) phase close to the transition point of the L = 3L and L = 5L surfaces, where 
iV = 1442, iV = 2562, and iV = 4842. The straight line drawn on the empty symbols was 
obtained by fitting the data to Eq. (JSJ) excluding the largest data, and the line indicates 
that a is constant independent of 2L/N in the smooth phase. Thus, we have z/~0 in 
the smooth phase. On the contrary, the exponent 

v = 0.93 ±0.14 (9) 

was obtained in the crumpled phase, and this indicates that the string tension is expected 
to vanish in the crumpled phase in the limit of t— >0. The vanishing a is the property 
that a decreases with decreasing t (increasing N under fixed L(N)/L Q (N)) or that a 
increases with increasing L under fixed N in the crumpled phase. Thus, the observations 
(i), (ii), (iii) on a mentioned above can be understood form the scaling property shown 
in FigEJ 

Table 1. The curvature coefficient a and the length of boundary 2L, where a in Fig[9] 
were obtained. 



N (L/L ) 


a (crumpled) 


a (smooth) 


2L 


1442 (3) 


2200 


2300 


42 


2562 (3) 


2000 


2200 


55.8 


4842 (3) 


1900 


1950 


76.8 


1442 (5) 


2500 


2600 


70 


2562 (5) 


2200 


2400 


93 


4842 (5) 


2000 


2100 


128 



Table [T] shows the curvature coefficients a in the crumpled phase and the smooth 
phase, and the boundary length 2L, where a in FigJH]were obtained. 2L = 42 in the first 
row of the table implies that L (N = 1442) = 7, which is identical to that described in 
Section [3J 

The first-order transition should be reflected in S3/N such that S3/N 
disco ntinuously changes at the transition point. In order to see this, we plot S3/N 
in Figsril3l(a)- fT07 c) . We find that 5*3 has a jump or a gap under the conditions 
L(N) = 2L (N), L(N) = 3L (N), and L(N) = 5L (N). The discontinuity of S 3 /N 
in the N = 4842 surface is more clear than those in the smaller surfaces, and this also 
shows that the model undergoes a first-order transition. 

We should note that a possibility of continuous transition is not completely 
eliminated, although Si/N, a, and S^/iV^ appear to change disco ntinuously in our 
simulation data. In order to see this more clearly, we need more simulation data at the 
transition point. Then, we can find the order of the transition precisely by using the 
finite-size scaling analysis. 
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Figure 10. The intrinsic curvature energy S3/N vs. a obtained under the conditions 
(a) L(N) = 2L (N), (b) L(N) = 3L {N), and (c) L(N) = 5L (N). Solid lines were 
drawn by the multihistogram reweighting technique. 



5. Summary and conclusions 

We have investigated an intrinsic curvature model on dynamically triangulated spherical 
surfaces by fixing two vertices as point boundaries separated by the distance 2L(N). 
Three different values are assumed for L(N); L(N) = 2L (N), L(N) = 3L (N), and 
L(N) = 5Lq(N), where Lq(N) is the radius of a sphere of size N and is chosen such that 
Si/N is almost identical to Si/N —1.5. The surface size are iV = 1442, iV = 2562, and 
iV = 4842. 

A deficit angle term S 3 is assumed for the curvature energy, and the Gaussian bond 
potential Si is included in the Hamiltonian; S = Si + a S3, where a is the curvature 
coefficient. The model is expected to have the smooth phase at sufficiently large a on 
dynamically triangulated surfaces. It is also expected that the smooth phase is separated 
from a non-smooth phase by the first-order transition, where the surface is almost linear 
and is very similar to the branched polymer surface [47]. We also confirmed in [43] 
that the order of the transition changes from first-order to second-order on tethered 
surfaces as the distance L(N) increases. Therefore, we aimed in this paper to see how 
the transition changes depending on L(N) on dynamically triangulated fluid surfaces, 
where the surfaces are expected to be oblong at sufficiently large L(N). 

It was confirmed in this paper that the first-order transition remains unchanged 
and is not influenced by how large L(N) is, up to L(N) = 5Lq(N) at least. In the 
crumpled phase, the surface spans from one boundary point to the other as an oblong 
and almost linear surface under the conditions 2Lq(N) <L(N) <5Lq(N). On the other 
hand, there appears a spherical lump on the surface between the two boundaries in 
the smooth phase. Almost all vertices are included in the lump, and therefore physical 
quantities, such as the bending energy S2/NB, change depending on whether the lump 
is included in the surface or not. A discontinuous change of 5*2, which is not included in 
the Hamiltonian, characterizes the first-order transition. The deficit angle term in the 
Hamiltonian also changes discontinuously at the transition point. 

The string tension a was obtained when L(N) = 3Lq(N) and L(N) =5Lq(N), where 
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the surface is sufficiently oblong even in the smooth phase though the lump is included 
in the surface. We found that a has a jump at the transition point, and that a takes a 
constant value in the smooth phase in the sense that it is independent of the distance 
L(N). On the contrary, a increases with increasing L(N) in the crumpled phase. 

The scaling property a oc (2L/N) U was also obtained at the smooth phase and at 
the crumpled phase close to the transition point. We have z/~0 in the smooth phase, 
and z/ = 0.93 ± 0.14 in the crumpled phase. Thus, the crumpled phase is characterized 
by vanishing string tension; a decreases with decreasing t = 2L/N (increasing N under 
fixed L(N)/Lq(N)) or equivalently a increases with increasing L under fixed N in the 
crumpled phase. This describes the above mentioned dependence of a on L(N) in the 
crumpled phase. 

It is interesting to calculate a in the extrinsic curvature model on relatively large 
fluid surfaces. The shape of surfaces of the extrinsic curvature model [28] seems 
very different from those of the intrinsic curvature model, and the dependence of the 
transition on L is also different from each other; the transition strengthens (remains 
unchanged) as L increases in the extrinsic (intrinsic) curvature model. Moreover, the 
dependence of a on the curvature coefficient in both models is almost identical with each 
other; a has a jump at the transition point under sufficiently large L in both models. 
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